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Abstract 

A systematic theoretical evaluation of the cross section for hadron 
production in electron-positron collisions in the energy range just be- 
low the B meson threshold is presented which includes charm and 
bottom mass effects and is accurate to order a^. The corresponding 
measurement in the energy region several GeV above the threshold is 
also discussed. 
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The hadronic decay rate of the Z boson provides one of the most accurate 
values for the strong coupling constant a s . This measurement is free from un- 
certainties and ambiguities which are inherent in the precise determination of 
a s from shape variables like thrust distribution, energy-energy correlation or jet 
multiplicities and which originate from the hadronization models, ambiguities in 
the jet definition etc. With increasing statistics and precision at LEP the uncer- 
tainty in a s from this measurement alone can be reduced to ±0.002. It would 
be highly desirable to test the evolution of the strong coupling as predicted by 
the beta function through a determination of a s from essentially the same ob- 
servable, however, at lower energy. The region from several GeV above charm 
threshold (corresponding to the maximal energy of BEPC around 5.0 GeV) to 
just below the B meson threshold at around 10.5 GeV corresponding to the "off 
resonance" measurements of CESR seems particularly suited for this purpose. As 
a consequence of the favorable error propagation, the accuracy in the measure- 
ment (compared to 91 GeV) may decrease by factor of about 3 or even 4 at 10 
and 5.6 GeV respectively, to achieve comparable precision in Aq C d- 



Most of the results derived in |l|, [| for massless quarks are applicable also for the 
case under consideration. However, two additional complications arise: 

(i) Charm quarks effects cannot be ignored completely and should be taken into 
consideration through an expansion in the ratio m 2 /s, employing the results of 
[0, |J for terms of order m 2 c /s and m\j s 2 . 

(ii) Contributions involving virtual bottom quarks are presented, starting from 
order a 2 . Their contribution depends in a nontrivial manner on m 2 /s. In order a 2 
these can be calculated in closed form and are shown to be small. Estimates for 
the corresponding contributions of order a\ indicate that these are under control 
and can be safely neglected, provided that one works within the correctly defined 
effective four quark theory. 

Through most of this paper the results will be formulated in the MS scheme 
for a theory with rif = 4 effective flavours and with the corresponding definitions 
of the coupling constant and the quark mass. The relation to a formulation with 
/ = 5 appropriate for the measurements above the bb threshold will be given at 
the end of the paper. 

We shall now list the independent contributions and their relative importance. 
Neglecting for the moment the masses of the charmed quark and a forteriori of 
the u, d and s quarks one predicts in order oP s 
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for the nonsinglet contribution [|I|, 0j. The second and the third order coefficients 
are evaluated with rif = 4 which means that the bottom quark loops are absent. 
In order a 2 the bottom quark loops can be taken into consideration with heir full 
mass dependence. However, as shown in || the leading term of order a 2 /ml, which 
has also been calculated in in the framework of an effective n/ = 4 theory, 
provides a fairly accurate description even up to the very threshold s 
Hence one has to add a correction term 
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For the singlet term one obtains 
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The bottom quark is absent in this sum. In view of the smallness of the a 2 s s/ml 
correction (even close to the bb threshold!) also all other terms of O(a^) from 
virtual b quarks are neglected. In the same spirit it is legitimate to use the scale 
invariant value of the b quark mass m& = m b (ml) as defined in the 5 quark theory. 
Starting from a pole mass of 4.7 GeV and using the results derived in || f| we 
get m b = 4.10 GeV if af\M z ) = 0.120. (Values of m b = 4.59 ± 0.04 GeV Jig, 
m b = 4.72 ± 0.05 GeV fill and m b = 4.7 ± 0.2 GeV |L2[] are currently in use for 
the pole mass. ) 

In contrast to the bottom mass the effects of the charmed mass can be in- 
corporated through an expansion in m 2 c / s. Quadratic mass corrections have been 
calculated up to order a z s1 quartic mass terms up to order a 2 s ||. Since m 2 J s is 
in itself a small expansion parameter, the order a 2 s m\j s 2 terms should be sufficient 
for the present purpose. 

The results for these corrections read 
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where rif — 4 has been adopted everywhere. Note that terms of order o? s m 2 c j s are 
more important than those of order a 2 m^/s 2 in the whole energy region under 
consideration. For completeness also m^/s 3 and a s m®/s 3 terms are listed, which, 
however, are insignificant amd will be ignored in the numerical analysis. 

The charm quark mass is to be taken as m c = frv£' (s) and is to be evaluated in 
the 4 flavour theory via the standard RG equation with the initial value m c (m c ) = 
1.11 GeV corresponding to a pole mass of 1.46 GeV in the case of a^(Mz) = 
0.120 0. 

A similar line of reasoning could be pursued for bottom mass terms in the 
region several GeV above the B meson threshold. It has been argued in j|, [13| 
( see also [Q) that mass effects of order a s are very well parametrized by taking 
leading terms of the expansion in m 2 q /s. This holds true not only in the high 
energy region but even a few GeV above the threshold. With this motivation in 
mind the terms of order a 2 ml/s and a s m^/s 3 will be included in the formula 
below. The result is expected to provide a reliable answer for y/s around 15 GeV 
and perhaps even down to 13 GeV. The corresponding formula reads as follows 
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Above the B meson threshold it is more conveniently to express all quantities for 
nf = 5 theory and thus in (^) all the coupling constant and quark masses are 
evaluated in the 5 flavour theory at the scale \i = y/s. 

The transition from the 4 to 5 flavour theory is performed as follows: The 
charm mass is naturally defined in the rif = 4 theory. In order to obtain the 
value of m c = rh^(s) the initial value m^lGeV) is evolved via the nj = 4 RG 
equation to the point /z 2 = m 2 and then from there up to fi 2 = s, however, now 
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with the rif — 5 RG equation. The bottom mass, on the other hand, is naturally 
defined in the rif = 5 theory irrespective of the characteristic momentum scale of 
the problem under consideration. As a numerical value we take the fh b {s) obtained 
from the scale invariant mass fhb(rhb) after running the latter with the help of 
the nf — 5 RG equation. 

Finally, and are related through the following equation 
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with x = ln(m^(/i)//i 2 ). Given Aq CD , eq. @ is employed to find A.q CD by setting 
ix = rn b (m b ). 

In tables 1-6 the predictions for R are listed for different values of af\Mz) 
together with the values of a s {s) and the running masses. (Note that our pre- 
dictions are presented without QED corrections from the running a and initial 
state radiation.) Figure 1 shows the behavior of the ratio R{s) as a function of 
energy below and above the bottom threshold, for a s (Mz) = 0.120,0.125 and 
0.130. The light quark (u, d, s, c) contribution is displayed separately also above 
10.5 GeV. It is evident that the predictions from the 4 and 5 flavour theories join 
smoothly. The additional contribution from the bb channel is presented down to 
11.5 GeV, where resonances start to contribute and the perturbative treatment 
necessarily ceases to apply. Evidently the bb channel is present with full strength 
down to the resonance region - an important consequence of QCD corrections. 
From our discussion it should be evident that the theoretical prediction is well 
under control. Mass effects are small below the bb threshold as well as a few GeV 
above. An experimental test is of prime importance. 
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Table 1: Values of AgL, AgL «( 4) (s), m^{s) and m b (m b ) &t y/s = 5.0 GeV for 
different values of af\M"l) ■ 



a s \ 1V1 Z) 


a(5) 
A MS 


A(4) 
A MS 




m c [S) 
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zoo ivie v 


Q1 7 MpV 

o 1 1 ivie v 


U.Zlo 


n 79^ PoV 
U. 1 ZO Vore V 


4 1 n PpV 


0.1225 


266 MeV 


357 MeV 


0.227 


0.666 GeV 


4.07 GeV 


0.1250 


302 MeV 


399 MeV 


0.236 


0.599 GeV 


4.04 GeV 


0.1275 


341 MeV 


444 MeV 


0.245 


0.522 GeV 


4.01 GeV 


0.1300 


383 MeV 


493 MeV 


0.255 


0.433 GeV 


3.98 GeV 



Table 2: Predictions for R(s) at = 5.0 GeV; the contributions to 5R mc are 
shown separately for every power of the quark mass. 





Rns 


Rs 


5R m 2 


5R m i 


SRm b 


R 


0.1150 


3.558 


-0.00015 


0.05 


-0.0066 


0.0012 


3.603 


0.1175 


3.567 


-0.00016 


0.047 


-0.0052 


0.0013 


3.610 


0.1200 


3.576 


-0.00018 


0.044 


-0.0040 


0.0014 


3.618 


0.1225 


3.586 


-0.00021 


0.040 


-0.0028 


0.0015 


3.624 


0.1250 


3.596 


-0.00023 


0.034 


-0.0018 


0.0017 


3.630 


0.1275 


3.606 


-0.00026 


0.028 


-0.0011 


0.0019 


3.634 


0.1300 


3.617 


-0.00029 


0.020 


-0.0005 


0.0020 


3.638 
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Table 3: Values of AgL, AgL, «i 4) (s), m[ 4 )(s) and m b (m b ) &t y/s = 10.5 GeV 
for different values of af\M"l) . 



a s \ m Z) 


a(5) 
A MS 


A(4) 
A MS 


a s {S) 






U. 1 iou 


1 7^ MpV 
i i o ivie v 


OAR MpV 

z^tu ivie v 


U. 1UU 


n 7^Q PpV 
u. i oy Vore v 


A 1 fi PpV 
1U *ore V 
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zoo ivie v 


Q1 7 MpV 

o 1 1 ivie v 


n 1 7fi 

U. 1 ( 


n fiAA PpV 


4 i n PoV 


0.1225 


266 MeV 


357 MeV 


0.182 


0.589 GeV 


4.07 GeV 


0.1250 


302 MeV 


399 MeV 


0.187 


0.527 GeV 


4.04 GeV 


0.1275 


341 MeV 


444 MeV 


0.193 


0.457 GeV 


4.01 GeV 


0.1300 


383 MeV 


493 MeV 


0.199 


0.376 GeV 


3.98 GeV 



Table 4: Predictions for R(s) at = 10.5 GeV; the contributions to 5R rric are 
shown separately for every power of the quark mass. 



af){Ml) 


Rns 


Rs 


5R m 2 


5R m 4 


8Rm b 


R 


0.1150 


3.518 


-0.000081 


0.0068 


-0.00022 


0.0022 


3.526 


0.1175 


3.523 


-0.000089 


0.0063 


-0.00017 


0.0024 


3.532 


0.1200 


3.530 


-0.000097 


0.0057 


-0.00013 


0.0026 


3.538 


0.1225 


3.536 


-0.00011 


0.0050 


-0.000088 


0.0028 


3.543 


0.1250 


3.542 


-0.00012 


0.0042 


-0.000056 


0.0030 


3.549 


0.1275 


3.548 


-0.00013 


0.0033 


-0.000032 


0.0032 


3.555 


0.1300 


3.555 


-0.00014 


0.0023 


-0.000015 


0.0034 


3.560 
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Table 5: Values of AgL, A^,af\s), m^(s) and mf\s) &t y/s = 13 GeV for 
different values of aj, 5) (M§ ) . 
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(5)/ \ 
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1 7^ MoV 
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3 47 HpV 
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u. izuu 
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n i 79 
u.± / z 


U.OZO Vore V 


^ 49 PoV 
0.41Z (ore V 


0.1225 


266 MeV 


0.177 


0.572 GeV 


3.36 GeV 


0.1250 


302 MeV 


0.183 


0.511 GeV 


3.30 GeV 


0.1275 


341 MeV 


0.188 


0.443 GeV 


3.23 GeV 


0.1300 


383 MeV 


0.194 


0.364 GeV 


3.17 GeV 



Table 6: Predictions for R(s) at =13 GeV; the contributions to 5R m are 
shown separately for every power of the quark masses. 



af){Ml) 


Rns 


Rs 


SR m 2 


SR m 4 


SR m 6 


R 


0.1150 


3.863 


-0.000019 


0.027 


-0.012 


-0.0016 


3.876 


0.1175 


3.869 


-0.000021 


0.027 


-0.011 


-0.0015 


3.884 


0.1200 


3.875 


-0.000023 


0.027 


-0.011 


-0.0014 


3.891 


0.1225 


3.882 


-0.000025 


0.027 


-0.0099 


-0.0013 


3.898 


0.1250 


3.888 


-0.000027 


0.027 


-0.0092 


-0.0011 


3.905 


0.1275 


3.895 


-0.000030 


0.027 


-0.0086 


-0.0010 


3.912 


0.1300 


3.902 


-0.000032 


0.026 


-0.0079 


-0.00091 


3.919 
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Figure 1: The ratio R{s) below and above the b quark production threshold at 
10.5 GeV for a s (M z ) = 0.120, 0.125 and 0.130. The contributions from light 
quarks are displayed separately. 
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